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We investigate localization properties of electron eigenstates in one-dimensional (Id) systems 
with long-range correlated diagonal disorder. Numerical studies on the localization length ^ of 
eigenstates demonstrate the existence of the localization-delocalization transition in Id systems and 
elucidate non-trivial behavior of ^ as a function of the disorder strength. The critical exponent 
u for localization length is extracted for various values of parameters characterizing the disorder, 
revealing that every v disobeys the Harris criterion v > 2/d. 
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I. INTRODUCTION 



Spatial correlation of disorder often causes an un- 
expected phenomenon in quantum disordered systems. 
Most intriguing is the breakdown of Anderson localiza- 
tion in one-dimensional (Id) systems induced by corre- 
lated disorder. Breakdown of the localization in Id sys- 
tems has beeiipredicted for the first time for a random 
dimer model |]|, wherein the on-site potential {e;} has 
a binary distribution with a short-range spatial correla- 
tion. Subsequently, a discrete number of extended eigen- 
states was found numerically in the random dimer model 
Q, Q . It was examined experimentally using transmis- 
sion measurements on semiconductor superlattices 3. 
These findings have motivated the studies of the na- 
ture of Id systems with long-range correlated disorder 
UlllllillllSIIllil- Particularly noteworthy is the 
system in which the sequence {e^} has a power-law spec- 
tral density of the form S{k) (x k~^: S{k) is the Fourier 
transform of the spatial correlation function (eiSj). For 
exponents p greater than 2.0, there is a finite range of 
energy values with extended eigenstates [||. This fact 
indicates the presence of the localization-delocalization 
transition in Id disordered systems a gain st the conclusion 
of the well-known scaling theory 0, hj. The emergence 
of extended eigenstates was also observed in harmonic 
chains with random couplings and in those with ran- 
doms masses jl^J. Note that diagonal disorder treated 
in Ref. is characterized by the power-law spectral 
density denoted above. 

The randomness of the long-range correlated poten- 
tials {si} is characterized by two quantities. The first is 
the exponent of the power-law spectral density, p, deter- 
mining the roughness of potential landscapes, as shown 
in Fig. 1. The second quantity is the distribution width 
W defined by the relation S [—W/2, W/2] ^ which char- 
acterizes the amplitude of the potential [13 . Effects of 
p on the localization properties of eigenstates have been 
examined but those of W remain unclarified. If the 
disorder is spatially uncorrelated (p = 0), an increase in 
W trivially induces strongly-localized wavefunctions be- 
cause all eigenstates are localized exponentially 0, . 



On the other hand, when the disorder is sufficiently long- 
range correlated to yield extended eigenstates {p > 2.0), 
the system shows a critical point Wc separating localized 
and delocalized phases. However, there is no attempt 
to quantitatively determine the value of Wc- Calcula- 
tions for Wc for various values of p allow us to establish 
the phase diagram in the W-p space, thereby engender- 
ing better understanding of localization properties of the 
system. We are also interested in the value of the criti- 
cal exponent v for the localization length of eigenstates. 
Those values can be obtained accurately using finite-size 
scaling analysis. 

The present work is intended to reveal critical proper- 
ties of electron eigenstates in Id systems with long-range 
correlated disorder. Numerical studies on the localization 
lengths have demonstrated the non-trivial behavior of 
^ as a function of W and p. A series of critical widths Wc 
and that of critical exponents v are determined by finite- 
size scaling analysis. Remarka bly, the results of v disobey 
the Harris criterion, v > 2/d llj,!!^, which is believed 
to be satisfied in general disordered systems. Our find- 
ings introduce new prospects for the study of Anderson 
transition in Id systems with correlated disorder. 

This paper is organized as follows. Section II describes 
the long-range correlated on-site potentials to be consid- 
ered. It presents a numerical algorithm for calculating 
the localization length of eigenstates. Section HI ana- 
lyzes the localization length as a function of W and p. 
The transition point Wc and the critical exponent v are 
extracted by finite-size scaling analysis. Conclusions and 
discussion are presented in Sec. IV. 



II. MODEL AND METHOD 

A. Long-range correlated potentials 

We consider noninteracting electrons in Id disordered 
systems within a tight binding approximation. The 
Schrodinger equation of the system is expressed as 



(1) 
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expressed by the following matrix equation: 



2000 4000 6000 8000 

Site index 

FIG. 1: Landscapes of spatially correlated on-site potentials. 
The exponent p of the power- law spectral density S{k) <x. 
is varied as displayed in the figure. The roughness of the 
energy landscape is gradually reduced with increasing p. 



where (pi is the amplitude of the wavefunction at the «-th 
site of the lattice. The hopping energy t is taken as a unit 
of energy hereafter. A sequence of long-range correlated 
potential {£») is produced by the Fourier filtering method 
|23ll23 . l25| . This method is based on a transformation of 
the Fourier components of a random number sequence. 
The outline of the method is: i) A sequence {ui] of un- 
correlated random numbers with a Gaussian distribution 
is prepared, ii) Its Fourier components {uq} are com- 
puted using the fast-Fourier transformation method, iii) 
A sequence [Sq] is newly generated for a given p using 



the relation £„ 



-p/2, 



iv) Finally, the objective {si} 
is obtained as the inverse Fourier transform of {Sq}- The 
resulting potentials Si are spatially correlated and pro- 
duce the power-law spectral density S{k) oc . In the 
following, the mean value (e;) is set to be zero and the 
periodic boundary condition is imposed. 

Figure 1 displays the landscapes of long-range corre- 
lated potentials Si that are generated by the procedure 
above. The system size L — 2^^ and the distribution 
width W = 3.0 are fixed. An increase in p markedly 
reduces the roughness in potential landscapes. We have 
confirmed that the sequences {si] appearing in Fig. 1 
produce the power-law spectral density, S{k) cx fc"^, for 
all values of p. 



B. Localization lengths 



Localization lengths of eigenstates in a potential field 
{si} are computed easily using the conventional trans- 
fer matrix method (26l|. The Schrodinger equation is 



M, 



E-e, -1 
1 



(2) 

The localization length ^ at a given energy E is defined 
by the relation |2^ 



1 |n,^,M.z(o)| 



with a generic initial condition 



z(0) 



(3) 



(4) 



Equation gives the single value of ^ only for the infi- 
nite system size L — > oo. However, when L is finite, the 
calculated result of Eq. (jj)) - denoted as - depends on 
the choice of the potential field {si}. To obtain a typical 
value of for a given L, we take a geometrical mean of 
on more than 10'* samples. Energy E is fixed at the 
band center E — Q throughout this paper. 

A critical point W = Wc (and p — pc) can be de- 
duced from the dependence of the normalized localiza- 
tion length A = ^i/L on the system size L. The typical 
values of increase with L for delocalized states, where 
the growth of is faster than that in L. This causes 
the quantity A to be an increasing function of L. On 
the other hand, A vanishes for sufficiently large L for lo- 
calized states, because approaches a constant value. 
Therefore, at the localization-delocalization transition, A 
must be invariant for the change in the system size L. 
Values of Wc are obtained accurately using the finite-size 
scaling method, as explained in the next section. 



III. NUMERICAL RESULTS 
A. The M^-dependence for the function A{W) 

Figure 2 plots the normalized localization length A as 
a function of the distribution width W. The exponent p 
of the power-law spectral density S{k) is increased from 
p = 1.0 to p = 4.5 incrementally. The system size L = 2^^ 
is fixed for all p. When p equals unity or less, A(iy) 
is a monotonous function of W. For larger p, on the 
other hand, curves of A show a kink at W = 4.0, which 
sharpens as p increases. 

We find two striking features in Fig. 2. The first is 
a peculiar p-dependence of A{W). The p-dependence of 
A(M^) in the region where W < 4.0 differs completely 
from that in the region where W > 4.0. For W < 4.0, 
the values of A rise with increasing p indicating that the 
growth of p ioT W < 4.0 causes an increase in the local- 
ization length of eigenstates. In contrast, the values 
of A ioi W > 4.0 systematically decrease with increasing 
p (except for the data p = 1.0). Hence, growth of p for 
W > 4.0 produces strongly-localized eigenstates. This 




w 

FIG. 2: The normalized localization length A = ^l/I/ as a 
function of the distribution width W. Several values of p are 
taken as denoted in the figure. The system size L is fixed to 
be 2^'' for all values of p. 
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FIG. 3: The L-dependence for A{W) with setting p = 1.0 
(open symbols) and p — 2.5 (solid symbols). Inset: Enlarge- 
ment of A(W^) in the vicinity of VK = 4.0. The system size 
L is varied to L = 2^^ (solid circles), 2^* (open circles), 2^^ 
(solid triangles), 2^^ (open triangles), and 2^^ (diamonds). 



difference in the effect of increasing p on is non-trivial 
because the increase in p simply smoothes potential land- 
scapes, as shown in Fig. 1. Therefore our results suggest 
that the effect of potential roughness on the localization 
length depends strongly on the value of W. 

The second notable feature is a shoulder structure of 
A{W) - a sharp bent of A{W) ~ at around W = 4.0, 
engendering a plateau-like shape within the region 0.5 < 
W < 4.0. The shoulder of A{W) appears for exponents 
p > 2.0, i.e., for p large enough to yield extended eigen- 
states That fact implies that localization properties 
of eigenstates for W < 4.0 (and p > 2.0) differ substan- 
tially from those for W > 4.0. The value W = 4.0 is a 
critical disorder strength separating a localized and de- 
localized phase, as explained later. The next subsection 
presents a demonstration that eigenstates are delocalizcd 
within the plateau region, under conditions W < 4.0 and 
p > 2.0. 

Values of the critical width Wc can be estimated from 
the dependence of A{W) on system size L. Figure 3 
shows the L-dependence for A{W), where the system 
size is increased from L ~ 2^^ (circles) up to 2^^ (dia- 
monds). Open and solid symbols correspond to the ex- 
ponent p = 1.0 and p = 2.5, respectively. For p = 1.0, 
the magnitude of A(W^) declines with increasing L such 
that the eigenstates are localized for any W. On the 
other hand, for p ~ 2.5, the i-dependence for A(W^) is 
rather complicated. The inset of Fig. 3 displays detailed 
behavior of A in the vicinity of = 4.0, where the sys- 
tem size is varied from L — 2^^ to 2^^ as denoted in the 
figure caption. All curves intersect on a single point at 
W = 4.0, suggesting the presence of the Anderson tran- 



sition atW = Wc = 4.0 for p = 2.5. 

We have confirmed that, when p > 2.0, A(VF) shows 
the same L-dependence as presented in the inset of Fig. 3. 
Intriguingly, all data of A{W) for p > 2.0 provide an 
identical value of Wc — 4.0. Therefore, we conclude that 
the critical distribution width Wc — 4.0 is independent 
oi p whenever p > 2.0. 

B. Finite-size scaling analysis 

Finite-size scaling analysis allows the determination of 
critical properties of the transition for L — > oo from data 
for finite L [2^ . This method stems from the hypothesis 
that the normalized localization length A close to the 
transition obeys the scaling law expressed by 

.nA^/(A), ,5) 

where is the localization length of eigenstates in an 
infinite system. The argument L/S^ao becomes much 
smaller than unity because ^oo diverges with obeying 
the form \W — Wc\~^ near the transition point 

W = Wc- This allows expansion of the scaling function 
as 

In A = aQ+ai\W-Wc\L^^'' + - ■ ■+an\W -Wc\'' L''^" , (6) 

terminating the expansion at the order n. Fitting the 
numerical data of In A for various values of W and L to 
Eq. (|SJ), we obtain Wc and v with great accuracy. Note 
that the optimal value of constants ag, oi, a„ for 
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FIG. 4: Scaling plots of In A = f{x) for p — 2.5 (solid circles) 
and for p — 3.0 (open circles). Upper and lower branches 
correspond to extended and localized phases, respectively. 



FIG. 5: The L-dependence for A(p) with setting W = 1.0 
(solid symbols) and W = 3.0 (open symbols), respectively. 
The system size L is varied as 2^'^ (circle), 2^^ (triangle), and 
2^"^ (rectangular). 



W > Wc are different from those for W < Wc because 
^oo is a function of the absolute value of Vt^ — Wc- 

Figure 4 shows scaling plots of hi A ior p — 2.5 (solid) 
and p = 3.0 (open). Here we define x = \W — Wc\L^^" 
and set n — 4. Each upper and lower branch in the figure 
corresponds to the extended phase (W < Wc) and the lo- 
calized phase {W > Wc), respectively. All data of In A 
for various values of W and L fit well onto two branches. 
Resulting values of the critical exponent v and the criti- 
cal distribution width Wc are listed in Table 1, where the 
error is a 95% confidence interval. We see that all values 
of Wc are almost identical to 4.0, as expected. On the 
other hand, values of v exhibit a discrepancy for different 
branches and different ps |23| . This discrepancy of v con- 
tradicts the principle of one-parameter scaling requiring 
that V should be independent of a choice of parameters 
in the Hamiltonian of the system. It is surmised that the 
discrepancy of v occurs because of a finite-sized effect 
that causes a systematic error in scaling plots. A novel 
tec hniq ue of scaling correction (two-parameter scaling) 
|28l 123 . l30| would help to solve the problem, as we shall 
present in a future study [SUI- 



TABLE L Calculated results of the critical exponent v and the 
critical distribution width Wc- The error is a 95% confidence 
interval. 
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C. The p-dependence for the function A(p) 



Next, we examine A behavior as a function of p. Figure 
5 shows a plot of the function A(p) with setting W — l.Q 
and 3.0, where the system size is varied from L = 2^^ 
to 2^^. For all W and L, A(p) increases with p\ it then 
becomes constant for p > 2.5. This p-dependence for A 
is consistent with the results displayed in Fig. 2, in which 
the growth of p for W < 4.0 causes an increase in A. De- 
tail calculations of A{p) in the vicinity oi p ~ 2.0 reveal 
that all curves of A(p) that belong to different Ls intersect 
at the point at p = 2.0 as long as W < 4.0. This indi- 
cates that localization-delocalization transition occurs at 
an identical point p = Pc = 2.0 whenever W < 4.0. The 
critical point pc = 2.0 estimated above is consistent with 
the conclusion reported in Ref. @. 

A series of critical points Wc and Pc we have found are 
summarized in the phase diagram illustrated in Fig. 6. 
Solid circles express the critical point W = Wc deduced 
from the finite-size scaling procedure, whereas solid rect- 
angulars express p — Pc defined by the position at which 
A{p) is independent of L (See Fig. 5). An extended phase 
appears in the region surrounded by the two straight lines 
of = 4 and p — 2. It is noteworthy that values of 
Wc and Pc equal an integral number: 4 and 2, respec- 
tively. The same integral value of Wc = 4 has been 
observed in other Id systems with correlated disorder 
[in Sin as discussed below. 
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FIG. 6: A phase diagram of the system in the W-p space. 
The critical line separating the localized and delocalized phase 
consists of two straight lines: W — Wc = 4.0 and p — Pc = 
2.0. 



IV. CONCLUDING REMARKS 

Two open problems remain with respect to the criti- 
cal properties of the transition in Id systems with long- 
range correlated disorder. First, all results of the critical 
exponent v listed in Table 1 disobey the Harris criterion 
V > 2/d 0, |23|. The inequality is widely believed to 
be satisfied in general disordered systems with any spa- 
tial dimension d, thereby determining the lower bound 
for V. According to the inequality, v in Id systems must 
be larger than 2, which disagrees with our results. We 
note here that the Harris criterion was originally derived 
for a spatially uncorrelated disorder ,21,] . Therefore, the 
relation v > 2/d may be violated in systems with long- 
range correlated disorder. In fact, for c/assica^ percolation 
model, the inequality must be modified in the presence 
of long-range correlation in the site or bond occupations 
[S^ill^l- To elucidate the lower bound of v for Id systems 
considered, we should generalize the argument in Ref. ^ 
for the correlated disorder producing the power-law spec- 
tral density S{k) oc k~P. 



Secondly, the critical disorder width Wc is exactly 
equal to an integral number Wc — 4. Furthermore, the 
value of Wc is independent of the potential roughness 
characterized by p. In fact, the same value of Wc = 4 
has been found in other correlated potential models. One 
example is the random dimer model Q, |^ |^ , where the 
site energies Ea and ef, are assigned in pairs. The dis- 
order width W of the system is defined by the relation 
Ea - Eb e [-W/2,W/2\. If W is less than the critical 
value Wc ~ 4, the localization length diverges obeying 
the form ^oo{E) oc E~'^, yielding the delocalized eigen- 
state at S = 0. Another example is the Harper model 
[Uls^, in which site energies are described by a periodic 
function — {W/2) cos(27riw) with an irrational number 
uj. This model undergoes the localization-delocalization 
transition at W = A. That is, all eigenstates become crit- 
ical at = 4 [s^. We conjecture that the coincidence of 
the value of Wc in the three different models is evidence 
of an unexplored universality of Wc in Id systems with 
spatially correlated disorder [s^ . 

In conclusion, we have investigated critical properties 
of Id electron systems subject to long-range correlated 
disorder. The normalized localization length A of eigen- 
states shows a non-trivial behavior as a function of W 
and p, which indicates a peculiarity of Id systems ex- 
hibiting the Anderson transition. Detailed calculations 
for A{W,p) reveal that the transition point Wc = 4 is in- 
variant to the change of the potential roughness. More- 
over, the critical exponent v in the considered system is 
proven to violate the Harris criterion v > 2/d. We hope 
that our findings enlighten the study of quantum phase 
transition in manifold disordered systems. 
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